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1. Introduction 

Understanding the phase diagram of QCD in the temperature - chemical potential (T,/x) plane 
is crucial for many implications in astrophysics, cosmology and in the phenomenology of heavy-ion 
collisions. While the lattice formulation of QCD provided fruitful information at finite temperature, 
at non- vanishing chemical potential lattice simulations suffer from the sign problem. Only in recent 
years several different methods Q, ||| 0] have been devised to at least partly overcome 

this obstacle. 

In this work we present results from lattice calculations in QCD with dynamical light and 
strange quark degrees of freedom at non zero temperature and imaginary chemical potential. Our 
calculations are performed with a tree level Symanzik-improved gauge action and an improved 
staggered fermion action, the p4-action with 3-link smearing (p4fat3) - We focus on a range of 
temperatures in the vicinity of the pseudocritical temperature T c at vanishing chemical potential, 
0.94 < T/T c < 1.08. At each temperature we carried out simulations on lattices with temporal 
extent N t = 4. Following []9|], at each temperature the strange quark mass was adjusted to its 
physical value and the light up and down quark masses were taken to be degenerate and equal to 
m s /!0, which corresponds to a constant Goldstone pion mass of about 220 MeV. This allows us to 
utilize the zero temperature results and interpolations of [§] for the conversion of lattice parameters 
into physical units. For definiteness, we apply a value of T c = 202 MeV for the pseudocritical 
temperature at zero density. To insure small finite volume effects the spatial volume has been 
chosen to be V l ^T = 4. The number of gauge field configurations analyzed varies from 1000 
to 2000, the configurations are separated by 10 trajectories. All numerical simulations have been 



performed using the Rational Hybrid Monte Carlo (RHMC) algorithm [ |10j , |11| ]. 

At each temperature we performed calculations at several values of the imaginary quark chem- 
ical potential, \L q = i^li,q = u,d,s which was taken to be degenerate for all three flavors. In the 
following, we first present results relevant for an estimate of the pseudocritical line in the (r,/^/) 
plane and then compare our findings at imaginary chemical potential with results obtained within 
the Taylor expansion approach. 



2. The pseudocritical line at imaginary chemical-potential 

At imaginary chemical potential and high temperature, the phase diagram of QCD features the 
Roberge Weiss transition at }Xf = 7iT /3, associated with the phase of the Polyakov loop. The line of 
Roberge Weiss transitions ends at T = Trw At lower temperature the 'quark-gluon plasma' region 
is limited by a chiral transition [[| ||] which continues to real values of jl. In this work we focus 
our study on the range of temperatures 0.94 <T/T C < 1.08, i.e. to the vicinity of the pseudocritical 
temperature T c at zero density. 

In order to study the phase structure of QCD in the (T,Hj) plane, it is very useful [Q] to 
consider the phase of the Polyakov loop, L(x), that we can parameterize as L{x) = \L{x)\e"^ . In the 
presence of dynamical fermions and with imaginary chemical potential we expect (<p) = — 8 at low 
temperatures (6 = Hi/T), and (0) ~ 2kn/3 for (k— 1/2) < j^d < (k+ 1/2) at high temperatures; 
the values 8 = 2(^+1/2)^/3 correspond to the Roberge Weiss transitions from one Z3 sector to 
the other. In Fig. [j] we show our results for (0 ) versus the imaginary chemical potential for different 
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Figure 1: The Polyakov loop phase (<j>) as a function of the imaginary chemical potential for 4 different beta 
values. The vertical dashed line corresponds to the RW transition at a/ij = 7i/3N t . 



values of temperatures. For j8 = 3.290 which is below the critical value at jij = we found that 
(<j>) starts to deviate from zero right away at non vanishing potential. For the j8 values above the 
critical one we observe that (0) stays close to until a certain temperature dependent value of ajij 
is reached where the Polyakov loop phase begins to decrease with increasing jij. We will take that 
value as an estimate for the location of the pseudocritical chiral line [Q]. A jump of the phase is not 
observed at any temperature which indicates that all our T values are below Trw 

In Fig. || we display our results for the modulus of the Polyakov loop and the light quark chiral 
condensate as functions of the imaginary chemical potential for j8 = 3.320,3.335,3.351 which 
correspond to temperatures of 204.8 MeV, 209.6 MeV and 218.2 MeV respectively [§]. The chiral 
condensate is defined as 

(Y¥) q = \j^(TrM- l (m q )), q = u,d,s (2.1) 

where N% is the spatial volume and M the fermionic matrix. In the same figure we also display the 
behavior of the Polyakov loop susceptibility 

XL =Nl((L 2 )-(Lf) (2.2) 

as a function of for the same beta values. In the upper row of Fig. || we indicate by vertical lines 
the estimates for the critical /I/ values and their errors as obtained from the phase of the Polyakov 
loop. It is seen, most clearly at the largest temperature investigated, that in the region where the 
Polyakov phase starts to deviate from also the chiral condensate exhibits a rise while the modulus 
of the Polyakov loop decreases. 

Based on the estimates for the pseudocritical }Xj values, to leading order in fx/T the critical 
line can be parameterized as 




(2.3) 
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Figure 2: Chiral condensate and Polyakov loop as functions of the imaginary chemical potential (upper 
row) for beta values above the critical value. The vertical dashed line corresponds to the RW transition, 
the vertical continuous ones to the estimation of the critical jij from the Polyakov loop phase; (lower row) 
Polyakov loop susceptibility as function of fij for the same beta values. 



where T c is the critical temperature at zero chemical potential. The coefficient of the leading term, 
?2, has been calculated for various cases, see the collection in [12]. In Fig. ||] (left) we compare 
our data points on T c (jx)/T c with results which have been obtained from reweighting simulations 
at jit = within the same discretization scheme as ours. In [^] two rather heavy quark flavors 
were taken into account whereas in [14] the number of flavors was 3, with similar masses as in our 
case. Although the error bars are large the comparison suggests that the curvature t% grows with Nf 
which is consistent with a behavior ~ Nf/N c found in large N c expansion [15]. In Fig. || (right) we 
show our results fit. We checked the stability of the fit by choosing different ranges in the chemical 
potential values. Considering the entire range of values we find t% = 0.89(4), discarding the first 
value we find t% = 0.89(5) and discarding the last one, ti = 1.02(12). 
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Figure 3: The pseudocritical line at imaginary pi: comparison between our data points and results from the 
literature |§] (left); our fit results (right). 
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3. The quark number density at imaginary chemical potential 



For a large homogeneous system, the pressure and its first derivatives, the quark number den- 
sities, are defined as 



1 



lniT 



P_ 



VT 3 
1 dln3f 



(3.1) 
(3.2) 



where the partition function 3f is a function of the volume V, temperature T, quark masses m u ^, s 
and chemical potentials jJ. U:C i,s- Note that for imaginary chemical potential the quark number density 



is purely imaginary Qlq]. 

We calculated the quark number density at various values of the imaginary chemical potential 
/i = ijJ,j. In Fig. |] (left) we show our results for the light, «/, I = u,d, and the strange quark number 
density, n s , at the temperature T = 209.6 MeV. We fitted our data to the ansatz 



2\e 



3n(jU/) =A^/(c 2 -/l/) 



(3.3) 



which was suggested in [17]. This ansatz takes into account that the quark number density is 
an odd function of the chemical potential. Moreover, if the exponent e is less than 1, Eq. 3.3 



leads to a singularity in the quark number susceptibility at c. Eq. 3.3 can also be derived from 
the singular part of the ln2£ in the vicinity of a critical point c = with e given by the critical 
exponent a as e = 1 — a. In Fig. || (right) we show the results of the fit to the light quark density 
obtained at a temperature of 209.6 MeV. A fit to our entire interval and with c unconstrained gives 
A = 105(67), c = 0.285(14) and e = 1.30(28). In the range [0.15,0.26] of a\i u we obtain A = 
57(54), c = 0.276(15) and e = 1.06(36). We thus find that e is consistent with 1 and substantially 
larger than the value 0.28(2) found in [ |l7| ] while c exceeds %T /3 slightly. We interpret these 
results as to indicate that the regular contributions to ln?¥ are dominating the number density at 



the temperatures investigated (see also [13]) 
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Figure 4: (left) «/ and n s in comparision; (right) quark number density n\ at T=209.6 MeV fitted to the 
form predicted by a simple critical behavior at imaginary chemical potential. The vertical line indicates the 
Roberge Weiss transition. 
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Figure 5: The light quark number density at 209.6 MeV calculated at imaginary chemical potentials (points) 
and from Taylor expansion (curves with error bands). The pink curve from [ |l8| ] takes into account coeffi- 
cients up to 6th order. The blue curve has been obtained from our own computations of the coefficients until 
fourth order which have less statistics. 



One of the approaches to QCD at finite density is based on the Taylor expansion of the pressure 
in terms of the quark chemical potentials around a generic jio = (Hufi,Hd,o,Hs,o)> 



P_ 



(A) = £ c U n{P>u-hp) k {P>d-p>d$) l {$s-M n (3-4) 



k.Ln 



with the abbreviation \i q = \i q jT and the coefficients 

Ckln k\l\n\dfi u k dix d l dfi s n \T*) ^ 
evaluated at /in- Correspondingly, the number density of e.g. the u quark is given by 

Sf(A) = L kc kln{fru ~ fai,o) k ~ l (fld ~ Arf,o)'(A* ~ Ai,<>)"- (3-6) 
1 k,ln 

Usually the coefficients are computed at pL q Q = but they could also be calculated at imaginary 
values for the chemical potentials. This is work in progress. 

In Fig. H we compare our data for the quark number density calculated at imaginary chemical 
potentials with predictions from the Taylor expansion at a temperature of 209.6 MeV and at the 
same lattice parameters. Up to pLj/T ~ 0.4 there is practically no difference between the predictions 
obtained from Taylor expansion up to fourth and 6th order [18]. At larger values of pLj the 6th order 
Taylor curve starts to deviate from the fourth order one and bends downwards, thus qualitatively 
describing this trend of the data correctly. Still, the errors arising from trucanting the Taylor series 
at sixth order become sizeable at Hi/T ~ 0.4 at this temperature. 



4. Summary and conclusions 

We have studied QCD with 2+ 1 flavors at nonzero temperature and nonzero chemical potential 
on a line of constant physics with the strange quark mass adjusted to its physical value and the pion 
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mass of about 220 MeV. The simulations have been carried out at imaginary chemical potentials. It 
turned out that the temperatures of up to 1.08T C at which the computations were performed are be- 
low the endpoint Trw of a line of Roberge Weiss transitions. Instead, we could identify a crossover 
line at which phase and modulus of the Polyakov loop as well as the chiral condensate change their 
behavior. The curvature of this line was estimated as ?2 = 0.89(2) which is not inconsistent with 
the literature. Furthemore.., we calculated quark number densities at imaginary and compare our 
data with results obtained via Taylor expansion method. Up to \iijT ~ 0.4 good agreement was 
observed. 
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